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Abstract 



The chaotic scattering theory is here extended to obtain escape-rate expressions 
for the transport coefficients appropriate for a simple classical fluid, or for a 
chemically reacting system. This theory allows various transport coefficients 
such as the coefficients of viscosity, thermal conductivity, etc., to be expressed in 
terms of the positive Lyapunov exponents and Kolmogorov- Sinai entropy of a set 
of phase space trajectories that take place on an appropriate fractal repeller. This 
work generalizes the previous results of Gaspard and Nicolis for the coefficient of 
diffusion of a particle moving in a fixed array of scatterers. 
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1 Introduction 



One of the remarkable results of recent studies of non-equilibrium processes taking place 
in large systems has been the establishment of new connections between non-equilibrium 
statistical mechanics, and dynamical system theory . Although much more work needs 
to be done before these connections are clear and well established, some of the new results 
establish a direct connection between quantities of interest for the statistical mechanics 
of irreversible processes, such as transport coefficients, and quantities of interest for the 
dynamical description of the same system, such as Lyapunov exponents, and Kolmogorov- 
Sinai entropies. Here we will restrict our discussion to classical systems for which a reasonable 
level of understanding has been achieved. 

One method for developing such connections, which has been developed by Hoover, Posch 
and coworkers ||, and by Evans, Cohen, Morris and coworkers 0], is based on a study of 
systems that are maintained in non-equilibrium steady states by means of specially con- 
structed external or internal "forces" and a thermostat that removes energy generated in the 
system by the special forces. This method leads to interesting connections between transport 
coefficients, such as the coefficients of diffusion or viscosity, and the Lyapunov exponents of 
the thermostatted, forced system. For the case of a thermostatted single particle moving (in 
two dimensions) in a periodic Lorentz gas under the action of an electric field, the connec- 
tion between the diffusion coefficient, or the electrical conductivity, and the (two) Lyapunov 
exponents was studied by Baranyi, Evans, and Cohen ||, and rigorous results establishing 
this connection were obtained by Chernov, Eyink, Lebowitz and Sinai ||. One feature of 
this approach is that a Hamiltonian description of the system is not possible, due to the 
introduction of the thermostat and the question naturally arises as to whether a connection 
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of the sort described is possible for purely Hamiltonian systems. 

This question was answered in the affirmative, for the case of a classical particle moving 
in an environment of fixed scatterers, by Gaspard, Nicolis, and coworkers |7|, ||. They were 
able to show that the diffusion coefficient of the moving particle could be obtained in terms 
of the escape rate for the particle from a bounded region of the scatterers. This escape rate, 
in turn, is obtained - following arguments of Kantz and Grassberger 0], Bohr and Rand [TO 



Tel and coworkers Grebogi, Ott, and Yorke JT2[ - in terms of the positive Lyapunov 
exponents and KS-entropy that characterize the set of orbits of the moving particle that are 
trapped forever in the bounded region occupied by the scatterers. This set of "trapped" 
orbits form a fractal set of trajectories of the moving particle, and is referred to as the 



"repeller" ||13||. The relation between the diffusion coefficient, -D, for the moving particle and 



the dynamical quantities on the repeller is given by |7|] 



D = lim f-' 



]T \(L) - h KS (L) 

.Ai>0 



Here, we suppose that the scatterers are confined to a slab of width L in one dimension, 
and of infinite length in perpendicular spatial directions; ^ A . >0 Aj(L) is the sum over all 
positive Lyapunov exponents for trajectories on the repeller, and /iks(-^) is the KS-entropy for 
trajectories on the repeller. We assume, usually without proof, that the system has ergodic 
properties such that the Aj do not depend on the particle trajectory on the repeller. We also 
note that for the open systems considered here (i. e. particles not on the repeller escape 
from the region containing the scatterers), the KS-entropy on the repeller does not equal the 
sum of the positive Lyapunov exponents, and the difference is of order L~ 2 for large L, if the 
diffusion coefficient exists. The escape-rate expression for the diffusion coefficient has been 



evaluated explicitly for a two-dimensional periodic Lorentz gas ||, for a two-dimensional 
multibaker transformation WM , and for a variety of one- dimensional maps WE . 



In the preceding discussion, the diffusion process is actually considered as a chaotic 
scattering occuring in a large but finite system. Several recent works have been devoted to 
chaotic scattering and have elucidated the role of the fractal repeller in this phenomenon 
IE| . Usually, chaotic scattering is envisaged on a small scatterer containing a few scattering 



centers [II], |1J], [L7|]. However, it was shown that chaotic scattering becomes controlled by 
diffusion when the scatterer becomes large enough. Remarkable relations then exist between 
the properties of scattering and those of diffusion, which have recently been described in 



some detail [ITS 



It is the purpose of this paper to extend the chaotic scattering theory and, in particular, 
the escape-rate formalism to include other transport coefficients, and to include a treatment 
of chemical reaction-rate coefficients. This is accomplished by showing that the basic ideas 
used to derive the expression for the diffusion coefficient of a moving particle can be easily 
extended to apply to other transport and reaction-rate coefficients as well. 

In the next Sees. 2 and 3, we present the derivation of the escape-rate expressions for 
transport and reaction rate coefficients. This derivation is based on the fact that all of these 
coefficients can be related to the average mean square displacement of some appropriate 
dynamical quantity. In Sec. 4 of the paper, we outline some outstanding and interesting 
problems related to this work. A careful discussion of the theory for chemical reaction rates 
is given in Appendix A. Derivations of the relation between the escape rate of a dynamical 
quantity from a bounded region and the Lyapunov exponents and KS-entropy has been given 
in the literature for systems of a small number of dimensions |IJ, |II| [IT], |I2[ . In a separate 



paper JT9[, we will review this derivation, so as to extend it to systems with many degrees of 
freedom, to systems that may be of nonhyperbolic type, and to illustrate the method with 
applications to lattice gas automata and to hard sphere gas systems. 

2 Transport Coefficients and their Helfand Moments 

We begin by constructing a large system of N particles governed by Hamilton's equation of 
classical mechanics 

dH . dH , . 

« = W p = ~ W (2 ' 1) 

where (q, p) are the positions and momenta of the particles. We assume that these particles 
are contained in a rectangular domain of volume V. At the borders, we may consider either 
hard walls of infinite mass or periodic boundary conditions. In the latter case, the total 
momentum is conserved in addition to the total energy. Moreover, to simplify the relation 
with ergodic theory, we use the microcanonical ensemble and work on energy shells H = E. 

In the large system limit (N,V — > oo with N/V = n), irreversible processes in such 
a classical many-body system may be described by hydrodynamic equations such as the 
Navier-Stokes equations, or the diffusion equation, the chemical kinetic equations ruling the 
time evolution of chemical concentrations, or the equations of electrical conductivity which 
incorporate Ohm's law. These phenomenological equations contain dissipative terms which 
are dependent on transport and rate coefficients. The aim of nonequilibrium statistical 
mechanics is to obtain these coefficients in terms of the microscopic Hamiltonian equations. 
Since the work of Maxwell and Boltzmann, several methods have been developed to calculate 



these coefficients. The most general method, developed by Green and Kubo [£0], is the 



time correlation function approach whereby the transport and rate coefficients, are given as 
time-integrals of autocorrelations of the fluxes, 



a 



r lim (4 a) 4 a) ) dt , (2.2) 

Jo V^oo 

where is the flux at time t corresponding to the coefficient a. It is a function of the 
canonical variables (q,p) and is obtained by solving the equations of motion for a time t 
after an initial time so that 

4 a) (q,p) = 4 a) mq,P)} , (2-3) 

where denotes the flow in phase space induced by Hamiltonian equations (2.1). We can 
express this in terms of an iV-particle streaming operator 

4 a) = exp(-tL) 4 a) , (2.4) 

with the Liouvillian operator given in terms of a Poisson bracket expression L = {H, ■}. The 
average (•) in equation (2.3) is taken over a microcanonical ensemble in our case. 

It will be useful for us to apply the formulation of the Green-Kubo expressions, as ob- 



tained by moments G t such that the fluxes J t are derivatives, as [22 



4°° = j t Gi a) , (2.5) 



An integration by parts shows that 



((Gt ~ Go) 2 ) 



/ / {JtJt") dt' dt" = 2 t f (l - 
Jo Jo Jo V 



- ) (J J T ) dr . (2.6) 



Accordingly, if the following condition is satisfied 



1 ft 

lim - / r (JnJ T ) dr = 



we obtain the equality 



;2.7) 



lim 

t— >oo 



((Gt - Go) 2 ) 



2t 



POO 

/ (J J T ) dr . 
Jo 



(2.8) 



As a consequence, the transport and rate coefficients are also given by [22 



a = lim I lim ( - G q) ] 2 ) 



(2.9) 



Table I. Helfand's moments. 



process 


moment 








self-diffusion 


G (D) = 








shear viscosity 


G(v) = 


1 

VVk B T 


2-ii=l 


%i Piy 


bulk viscosity (?/> = ( + |t?) 


gw = 


1 


s^N 


%i Pix 


heat conductivity 


gw = 


1 


1 %i (E{ 


charge conductivity 


= 


1 

VVfc B T 




cZ{ Xi 


chemical reaction rate 


G (r) = 


1 

VVk B T 


(iV« 


- (iv (r) )) 



In Table 1 we list moments appropriate for each transport or rate coefficient, where 
^ = Jm I ^i(^i) ^ s ^ e ener gy °f particle z and V^j is the potential energy of interaction 
between particles % and j. Appendix A describes in more detail the case of chemical reactions. 
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Eq. (2.9) shows that, in the case the transport or rate coefficients are well defined (i.e. 
are positive and finite), Helfand's moments undergo a diffusive type of motion along the axis 
of the moment G^ a \ Therefore, the moment may be considered as a random variable having 
a probability density p(g) obeying a diffusion-type equation, in an equilibrium ensemble 

dp d 2 p 

m = a df 2 ' ( 0) 

with the transport or rate coefficient a a diffusion coefficient. This equation is the Fokker- 
Planck equation governing the equilibrium fluctuations of the Helfand moments. Equiva- 
lently, one can recover Eq. (2.9) by supposing that the moments satisfy a Langevin stochastic 
differential equation 



G< a > = 4 a) , (2.11) 

where the flux is a white noise 



(4 a) ) = , (4 a) 4 a) ) = 2a 8(t) . (2.12) 

It is important to note that equations (2.10) or (2.12) are to be considered as a simple 
representation of the results of the time correlation function method for time scales which 
are much longer than the time necessary for the time correlation functions of the microscopic 
currents to decay to zero. We suppose that this approach applies to situations where, for ex- 
ample, long time tails in the correlation functions decay sufficiently rapidly for the transport 
coefficients as defined by Eq. (2.2) to exist. In this way we can still treat processes which 
are diffusive on long time scales, but have correlations on shorter time scales. Accordingly, 
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the existence of a flux autocorrelation function which differs from a delta distribution on 
short time scales is still perfectly compatible with the validity of the diffusion-type equation, 
(2.10), on long time scales. 

In this discussion we have tacitly assumed that we consider a physical system in the 
proper thermodynamic limit. However, for systems of finite volume, there is an upper limit 
on the times we can consider because the ranges of variation of the Helfand moments are 
bounded. For example, if the system consists of hard spheres placed in a cubical box of 
length L on a side, the positions of each particle vary only in the interval 

L L . 

- - < Xi, yi, Zi < + - , (2.13) 

while the momenta can only take values in the interval 

- V2mE < p ix , p iy , p iz < + V2mE , (2.14) 

where E = E tot = ^NU-qT is the total energy. Similarly, the energy of a particle can only lie 
in the range 

< Ei < ^ k B T N . (2.15) 
Because of these bounds, the moments are always of bounded variation in the interval 

\G (a) \ < C N 5(a) , (2.16) 

where C is some constant, and 5^ are positive exponents which are respectively 5^ = 
1/3, = 4/3, <JW = 4/3, = 11/6, 6® = 5/6, and 5« = 1/2, for three-dimensional 
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gases. Accordingly, the range of variation of the moments grow with the size of the system, 
when the density and other intensive variables are kept constant. Of course, most variations 
of the Helfand moments will be due to microscopic motions of the particles and thus will 
be much smaller than the bounds in Eq. (2.16). At any rate, as the system gets larger, we 
expect that the diffusive-like behavior of the moments, described by Eq. (2.10), will be valid 
over increasingly larger regions of variations of the G^ a \ 

After the discussion about the limited range over which the Helfand moments obey the 
diffusive-type equation (2.10) for finite systems, we may proceed. 

3 Escape-Rate Formalism 

Within the range of validity of (2.10), we may set up a problem of first passage for the 
moment corresponding to the transport or rate processes of coefficient a. We consider 
a statistical ensemble formed by copies of the system which we assume to be at equilibrium 
and microcaninical at total energy E (and eventually at fixed total momentum in the case of 
periodic boundary conditions: P to t = Z)£Li Pi)- F° r eac h copy, the motion of the Hamiltonian 
system (2) is integrated from the initial conditions and the Helfand moment is calculated 
along the trajectory. At each time, we count the number of copies, J\f( a \t), for which the 
moment is still in the following interval 

- f < Gl a) < + | , (3.1) 

where the size of the interval x is sufficiently large to be in the regime of diffusion of the 
moment but not too large with respect to the total variation interval of the moment allowed 
by the fmiteness of the system. We are here defining a problem of first passage which can be 
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solved using the eigenvalues and eigenfunctions of the Fokker-Planck equation (2.10) with 
the boundary conditions 



p(- X /2) = K+X/2) = . (3.2) 
The solution of this eigenvalue problem is well known to be 



p(g,t) = f) c n exp(- 7 ^) sinf — g) , with 7 « = a (-Y , (3.3) 

where the constants c n are fixed from the initial probability density p(g,0). The number of 
copies of the statistical ensemble which are still in the interval (3.1) is then given by 



Af^(t) = [ +X/2 P (g,t)dg. (3.4) 

At long times, the decay is dominated by the slowest decay mode corresponding to the 
smallest decay rate 7l which defines the escape rate of the moment out of the interval 
(3.1), 



7T 



2 



li a) = a [-) . (3.5) 

We are now in position to establish a relationship with the deterministic dynamics. 
The Hamiltonian classical motion of the many-body system is chaotic in many cases. This 
property has been proved by Sinai and coworkers for some simple hard sphere gas models 
PB] , 53]. Also strong numerical evidence exists which shows that half of the Lyapunov 
exponents are typically positive in systems of statistical mechanics like the Lennard- Jones gas 



at room temperatures Po] . We suppose that the decomposition of phase space into ergodic 
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components is understood and that, beside the decomposition on the known constants of 
motion, there is only a single main ergodic component. 

We consider the set of all the trajectories for which the Helfand moment remains forever 
within the interval (3.1). Because most of the trajectories are expected to exit this interval, 
the trapped trajectories must be exceptional and highly unstable forming a set of measure 
zero with respect to the microcanonical probability measure. Based on earlier work on 
diffusion in the Lorentz gas and in related models and on the basis that the trajectories 
are typically of saddle type in systems of statistical mechanics |F| [8], [L4|, [L7|, [18 



we assume 



that the set of trajectories is a fractal repeller. Indeed, this set is of vanishing measure but 
may still contain an uncountable infinity of periodic and nonperiodic trajectories. A set 
satisfying these conditions is necessarily a fractal pi]. Moreover, it is composed of unstable 
trajectories of saddle type so that it forms a repeller (of saddle type) in phase space. These 



properties can be proved for particular models like the multibaker area-preserving map [14 



as well as the array of disk scatterers composing the periodic Lorentz gas. Accordingly, it 
seems reasonable to assume that, for more general systems, such as a gas of hard spheres 
the set of trajectories for which the moments satisfy Eq. (3.1), form a fractal repeller, with 
properties to be described in the next paragraph. 

A fractal repeller is characterized by different quantities and, especially, by an escape rate 
which is the deterministic analogue of the escape rate obtained in the preceding first-passage 
problem. Moreover, in chaotic systems, the escape rate is related to the sum of positive 
Lyapunov exponents minus the Kolmogorov- Sinai entropy per unit time if these quantities 
are well defined and positive |IJ, [| [II], O, [T7J . These quantities are evaluated for the natural 
invariant probability measure whose support is the fractal repeller. For the natural invariant 
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measure, each cell of phase space has a weight which is inversely proportional to the local 
Lyapunov numbers (stretching factors). Therefore we have [| [FT], |TT|, [T2|, [TJ], |T7 



l[ a) = E H^ a) ) - h KS (^) , (3.6) 

A»>0 

where we denote by J-^> the fractal repeller formed by the trapped trajectories for which 
the Helfand moment remains forever in the interval (3.1). 

Combining the deterministic result (3.6) with the statistical result (3.5), we obtain the 
relationship 



a = lim ( — | lim 



.Ai>0 



(3.7) 



where the limit V — > oo denotes the thermodynamic limit to be taken before the limit 
X — ► oo which is internal to the system. 

With Table I, Eq. (3.7) shows how a general transport or rate coefficient can in principle 
be related to the Lyapunov exponents and the Kolmogorov- Sinai entropy of a fractal repeller. 
This fractal repeller is the phase space object corresponding to the escape process of the 
Helfand moment associated to the transport or rate coefficient. In this way, a connection is 
established between statistical and mechanical considerations in phase space. 

A remark is now in order about the magnitude of the quantities appearing in (3.7). The 
sum of positive Lyapunov exponents and the KS-entropy per unit time are very large, of the 



order of the number of particles times the inverse of a typical kinetic time scale [p6] . On the 
other hand, the escape rate, which is the difference between two such large numbers, has a 
much smaller magnitude given by the time scales characteristic of hydrodynamics. In this 
way, the kinetic and hydrodynamic levels are naturally connected with a formula like (3.7). 

14 



4 Conclusion 

We have shown that all of the transport coefficients for a simple fluid and that chemical 
reaction rate coefficients can be expressed in terms of an escape rate from an appropriate 
fractal repeller. This completes a line of argument initiated by previous work on the coeffi- 
cient of diffusion for a particle moving in a periodic Lorentz system 0, and for a multibaker 
map [ 14| . In a separate paper [|T9 , we will show how the escape-rate formula can be ap- 



plied to several classes of dynamical systems like the hard sphere gas and the lattice gas 
automata. We also describe there a large deviation formalism which allows us to extend the 
application of the escape rate formula (3.7) from hyperbolic to nonhyperbolic systems, in 
particular, using the Ruelle pressure function. It now remains to apply this formalism to 
a number of examples in order to understand the physical and mathematical consequences 
of this approach to transport in fluid systems. Applications of this formalism to certain 



types of one- dimensional diffusion problems will be presented in [fLq1 , and applications to 
Lorentz lattice gas cellular automata will be presented in Many further applications 
are possible. 

Many interesting problems remain open in the present context, in particular: 

a) To provide experimental evidence for the microscopic chaos at the basis of the present 
theory. Suggestions along this line has been discussed elsewhere |28[ . 



b) To provide a more rigorous mathematical derivation of the escape-rate formalism used 
here. 

c) To establish the connection between this formalism and those based upon period-orbit 



theory ]29[] and upon a study of eigendistributions of the Liouville operator, and Ruelle 
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resonances flbj, |3(| pl| . Already, such a connection has been discussed elsewhere JT8 
where it was shown for the multibaker model that the diffusive-like eigendistributions of 
the Liouville operator shares self-similar properties with the fractal repeller underlying 
diffusion. Moreover, both approaches make use of the leading eigenvalues (3.5) of the 
Liouville operator and of its dependence on the wavenumber of the diffusive eigenmode 
(k = tt/L). 

d) To extend this formalism to include transport in mixtures, and to a treatment of higher 
order, and nonlinear transport processes. 

e) To study the structure of the fractal repeller for interesting cases where, as a con- 
sequence of long time tail effects, transport processes are anomalous, such as two 



dimensional fluid systems [EJ2] 



f) To establish the connection between the escape-rate formalism for transport and the 
Gaussian thermostat formalism of Hoover and Posch ||, and of Evans, Morris, Co- 
hen, and coworkers j|, ||. The Gaussian thermostatted systems belong to a class of 
dynamical systems which are time reversible without being volume preserving. The 
characteristic quantities describing chaos in these thermostatted systems have been 
described in [|, j33fl . 
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Appendix A: The Case of Chemical Reactions 

A.l. Summary of thermodynamic results 

For a chemical reaction like 



E < x i <- E ^ ^ , (A.i) 

7=1 7=1 

the numbers of particles of the reactants and products are changing at each step of the 
reaction according to 

A* = AN, = = ANc = Aiv(r) 
where z/ 7 = v ~ — z/+ are the stoichiometric coefficients [34]. The degree of progress of the 



reaction can be measured in terms of the variation of the number characterizing the reac- 
tion. We also introduce the chemical concentrations, C 7 = N 1 /V . At the phenomenological 
level of thermodynamics, the velocity of the reaction is defined by 



a W = °1 = W} (A.3) 



where C 7 are the average chemical concentrations ||34j| . The dependence of the reaction 



velocity on the concentrations themselves is given by the mass action law 



w = k + 1[ c? - n c ? ■ ( a - 4 ) 



7=1 7=1 

The affinity of the reaction is defined by 



A=-X> 7 ^ = - k B T In f[ (%T , (A.5) 

7=1 7=1 V °T/ 
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where /x 7 = /x° + In C 7 is the chemical potential of the species 7 and it is known that the 
affinity vanishes at thermodynamic equilibrium: A eq = Near the thermodynamic 



equilibrium, both the reaction velocity w and the affinity A can be expanded in terms of the 
variations of the chemical concentrations around their equilibrium values, C 7 = Cy 1 + AC 7 . 
In the linear approximation, we obtain the equality 



w ~ w c + q -4= , with = k + n ^ , (A.6) 

fc B^ 7= 1 

which allows us to obtain the Onsager coefficient of this chemical process 



LA . , w!_ q , , „. 

w ~ , with L = — . (A.7) 



Yamamoto [BO] and, later, Zwanzig p^] have shown that the Onsager coefficient of the chem- 



ical reaction is given by the following integral of the autocorrelation function in a classical 

system 



" (iV«iV«) dt , (A.8) 



Vk-Q Jo 

where the corresponding flux is here proportional to the time derivative of the number iV^ 
of particles which is characteristic of the reaction. Accordingly, we obtain the result given 
in Table I for the reaction coefficient 



a = — = — , (A.9) 
A T y ' 



with (JVW) = iV ( 



r),eq 
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A. 2. Master equation approach 

The escape-rate formula can be obtained for chemical reactions by using the master 
equation approach f55 |. As an example, we consider the isomerization 



A <-> B 



(A.10) 



The numbers of particles A and B are the random variables of this process. The total 
number of particles is conserved, N A + Nb = N = N tot , so that the process is completely 
determined by the knowledge of the number of particles A. The evolution equation of the 
probability P(Na) that the system contains Na particles is 



d 
dt 



P(N A ) = k + {N A + 1) P(N A + 1) + k-{N-N A + l) P{N A -1) 



k + N A P(N A ) - k-{N - N A ) P(N A ) . 



(A.ll) 



Introducing the fraction < / = Na/N < 1 of particles A and defining a probability density 
which is function of this continuous concentration according to p(f) = p(N A /N) = P(N A ), 
we use Taylor expansions in the master equation (A.ll) to obtain the Fokker-Planck equation 



^ + -Up) = - 
dt dfV p ) df 



D(f) ^ 



(A. 12) 



where 



f = - k + f + mi-/) + 0(1/N), 



(A.13) 



is the macroscopic rate equation while 
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D(f) = ^ [Kf + Ml-/)] , (A.14) 

is a concentration-dependent diffusion coefficient. This Fokker-Planck equation determines 
the time evolution of the probability density of a stochastic process in which the amplitude 
of the noise is proportional to 1/ \^N. Let us consider these thermodynamic fluctuations 
around the equilibrium concentration 



^ = • (A ' 15) 

in the vicinity of which the Fokker-Planck equation becomes a diffusion equation 



with 



D{r) = nKT^) ■ (A - 17) 

On this ground, we can apply the argument of first passage for the Helfand moment associ- 
ated with the chemical reaction, which is given in Table I 

G " = OT {N * " w?) ■ (A - 18) 

The diffusion coefficient (A. 17) shows that the moment is a random variable of Brownian 
type in a small neighborhood of the thermodynamic equilibrium. Indeed, the moment is 
related to the fraction / = N A /N of particles A according to 
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G( " = (/ " n ' (A - 19) 

We are looking for the first time at which the moment reaches the boundaries of the interval 



_ X < G (r) < + X _ ^ A ^ 20 ^ 
This first-passage problem corresponds to the first escape of the fraction / out of the interval 



r - | < / < r + \ , (A.21) 

with 



* = ^ X ■ (A.22) 

We consider a statistical ensemble formed by copies of the system and we count the 
number of copies which remain within the interval (A. 20) at time t. Since the moment 
or the fraction / follows a diffusion process, the number of copies is decaying exponentially 



A/" (r) (t) ~ exp - 7 { r) t , (A.23) 



(r) 

at a rate 7J given by 



* D(D , (A.24) 

in terms of the diffusion coefficient (A. 17) of the thermodynamic fluctuations around equi- 
librium. We mention that this exponential decay is the slowest decay dominating the time 
evolution at long times after the faster decay modes have died out. We can replace the size 
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e (A. 22) of the escape interval of the variable / by the size x °f the escape interval of the 
corresponding Helfand moment and we get 



7 J r) ~ a , (A.25) 



with the rate coefficient 



a = Vk B T(k + + k_) = A ' (A - 26) 
so that we recover the rate coefficient (A. 7) and (A. 9) of the macroscopic theory since 



„ = g, with „ ? = k+ C? = . (A.27) 

In this way, we see the consistency of the first-passage problem applied to the thermodynamic 
fluctuations described by the master equation with the macroscopic theory as well as with 
the Green-Kubo formula. The relation to the fractal repeller of the deterministic dynamics 
is discussed in full generality in the main part. 
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